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Abstract 

Motivated by biological questions, we study configurations 
of equal-sized disks in the Euclidean plane that neither pack 
nor cover. Measuring the quality by the probability that a 
random point lies in exactly one disk, we show that the regu¬ 
lar hexagonal grid gives the maximum among lattice config¬ 
urations. 

Keywords. Packing and covering, disks, lattices, Voronoi domains, 
Delaunay triangulations. 

1 Introduction 

High-resolution microscopic observations of the DNA or¬ 
ganization inside the nucleus of a human cell support the 
Spherical Mega-base-pairs Chromatin Domain model [3, 9], 
It proposes that inside the chromosome territories in eukary¬ 
otic cells, DNA is compartmentalized in sequences of highly 
interacting segments of about the same length [4]. Each seg¬ 
ment consists of roughly a million base pairs and resembles a 
round ball. The balls are tightly arranged within a restricted 
space, tighter than a packing since they are not rigid, and less 
tight than a covering to allow for external access to the DNA 
needed for gene expression. 

Motivated by these biological findings, [7] considered 
configurations in which the overlap between the balls is lim¬ 
ited and the quality is measured by the density, which we 
define as the expected number of balls containing a random 
point. We introduce a new measure that favors configura¬ 
tions between packing and covering without explicit con¬ 
straints on the allowed overlap. Specifically, we measure a 
configuration by the probability a random point is contained 
in exactly one ball. Since empty space and overlap between 
disks are both discouraged, the optimum lies necessarily be¬ 
tween packing and covering. The interested reader can find 
references for traditional packing and covering in two and 
higher dimensions in [2, 5], In this paper, we restrict atten¬ 
tion to equal-sized disks in the plane whose centers form a 
lattice, leaving three and higher dimensions as well as non¬ 
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lattice configurations as open problems. Our main result is 
the following non-surprising fact. 

Theorem 1 (Main) Among all lattice configurations in R, 
the regular hexagonal grid in which each disk overlaps the 
six neighboring circles in 30° arcs maximizes the probability 
that a random point lies exactly in one disk. 

For obvious reasons, we call this the 12-hour clock config¬ 
uration. We prove its optimality in four sections: preparing 
the background in Section 2, proving an equilibrium condi¬ 
tion in Section 3, developing the main argument in Section 
4, and giving the technical details in Appendix A. 

2 Background 

In this section, we introduce notation for lattices, Voronoi 
domains, and Delaunay triangulations. 

Lattices. Depending on the context, we interpret an el¬ 
ement of R 2 as a point or a vector in the plane. Vectors 

a, b £ R 2 are linearly independent if aa + fib = 0 implies 
a = fi = 0. A lattice is defined by two linearly independent 
vectors, a,b £ R 2 , and consists of all integer combinations 
of these vectors: 

L(a,b) = {ia + jb \ i,j £ Z}. (1) 

Its fundamental domain is the parallelogram of points aa + 
fib with real numbers 0 < a,fi < 1. Writing ||a|j for 
the length of the vector and 7 for the angle between a and 

b, the area of the fundamental domain is det L(a , b) = 
||a|||| 6 || sin('T'). The same lattice is generated by different 
pairs of vectors, and we will see shortly that at least one of 
these pairs defines a non-obtuse triangle. We will be more 
specific about this condition shortly, as it is instrumental in 
our proof of the optimality of the regular hexagonal grid. 

Voronoi domain. Given a lattice L, the Voronoi domain 
of a point p £ Lis the set of points for which p is the closest: 

v(p) ={reR 2 \\x-p\\ < 11*- g||,Vg e L}. (2) 

It is a convex polygon that contains p in its interior. Any two 
Voronoi domains have disjoint interiors but may intersect in 
a shared edge or a shared vertex. The lattice looks the same 
from every one of its points, which implies that all Voronoi 
domains are translates of each other: V(p) = p + V(0). 
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Figure 1: A primitive Voronoi diagram on the left, and a non¬ 
primitive Voronoi diagram on the right. The Delaunay triangula¬ 
tions are superimposed. 

Similarly, central reflection through the origin preserves the 
lattice, which implies that V (0) is centrally symmetric. 

The Voronoi diagram of L is the collection of Voronoi do¬ 
mains of its points. It is primitive if the maximum number 
of Voronoi domains with non-empty common intersection is 
3. In this case, the Voronoi domain is a centrally symmet¬ 
ric hexagon; see Figure 1. In the non-primitive case, there 
are generators that enclose a right angle, and the Voronoi do¬ 
mains are rectangles. To the first order of approximation, the 
area of any sufficiently simple and sufficiently large subset of 
R 2 is the number of lattice points it contains times the area 
of V(0). Similarly, it is the number of lattice points times 
the area of the fundamental domain. It follows that the area 
of V (0) is equal to det L. 

Packing and covering. For g > 0, we write B{p,g) 
for the closed disk with center p and radius g. The pack¬ 
ing radius is the largest radius, tl, and the covering radius 
is the smallest radius, Rl, such that B{0,tl) C V(0) C 
f?(0, Rl)- The density of the configuration of disks with ra¬ 
dius g centered at the points of L is the area of a disk divided 
by the area of the Voronoi domain: 

We) = O) 

It is also the expected number of disks containing a random 
point in ]R 2 . The packing density is Sl^l), which is neces¬ 
sarily smaller than 1. It is maximized by the regular hexag¬ 
onal grid, H, for which we have Sh^h) = 0.906 .... The 
covering density is Sl{Rl), which is necessarily larger than 
1. It is minimized by the regular hexagonal grid for which we 
have Sh(Rh) = 1.209 .... More generally, it is known that 
H maximizes the density among all configuration of con¬ 
gruent disks whose interiors are pairwise disjoint [10], and 
it minimizes the density among all configurations that cover 
the entire plane [8]. Elegant proofs of both optimality results 
can be found in Fejes Toth [5]. 

Delaunay triangulations. Drawing a straight edge be¬ 
tween points p and q in L iff V (p) and V (q) intersect along 
a shared edge, we get the Delaunay triangulation of L. In 
the primitive case, the edges decompose the plane into trian¬ 
gles. Among the six triangles sharing 0 as a vertex, three are 


translates of each other and, going around 0, they alternate 
with their central reflections. It follows that all six triangles 
are congruent and, in particular, they have equally large cir- 
cumcircles that all pass through 0. Since their centers are 
vertices of the Voronoi domain of 0, we have the following 
result. 

Lemma 2 (Inscribed Voronoi Domain) The vertices of 
V (0) all lie on the circle bounding B(0 , Rl)- 

The discussion above proves the Inscribed Voronoi Domain 
Lemma in the primitive case. It is also true in the simpler, 
non-primitive case in which V(0) is a rectangle. Return¬ 
ing to the primitive case, we note that the two angles op¬ 
posite to a shared edge in the Delaunay triangulation add 
up to less than 180°. In a lattice, these two angles are the 
same and therefore both acute. The two types of triangles in 
the Delaunay triangulation of a lattice can be joined across a 
shared edge in three different ways. We can therefore make 
the same argument three times and conclude that all angles 
are less than 90°. A slightly weaker bound holds in the non¬ 
primitive case. 

Lemma 3 (Non-obtuse Generators) Every lattice L in R 2 
has vectors a, b £ R 2 with L = L{a, b) such that 

(i) in the primitive case 0, a, b are the vertices of an acute 
triangle, 

(ii) in the non-primitive case 0,a, b are the vertices of a 
non-obtuse triangle with a right angle at 0. 

Assuming a, b satisfy the Non-obtuse Generators Lemma, 
the triangle 0 ab has edges of length ||a||, ||6||, and ||c|| = 
||a, — 6|L In the non-primitive case (ii), ||a|| and ||6|| are the 
lengths of the sides of the rectangle V(0), and ||c|| is the 
length of a diagonal. We have ||c|| 2 = ||a|| 2 + ||6|| 2 , and 
therefore ||a|| < ||6|| < ||c||, possibly after swapping a and 
b. In the primitive case (i), we can choose a, b, and c = a —b 
such that ||a|| < ||6|| < ||c||. In this case, V(0) is a centrally 
symmetric hexagon with distances ||a||, ||6||, ||c|| between an¬ 
tipodal edge pairs. 

3 Equilibrium Configurations 

Given a lattice in R 2 , we are interested in the radius of the 
disks for which the probability that a random point lies in¬ 
side exactly one disk is maximized. Further maximizing this 
probability over all lattices, we get the main result of this 
paper. 

Partial disks. Fix a lattice L in R 2 . For a radius g > 0, 
consider the set of points that belong to the disk centered at 
the origin but not to any other disk centered at a point of L: 

D(0,e) = B(0,g)\ (J B(p,e). 

ofpeL 


( 4 ) 
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As illustrated in Figure 2, for radii strictly between the 
packing radius and the covering radius, this set is partially 
closed and partially open. We distinguish between the con¬ 
vex boundary that belongs to the circle bounding B(0,g), 
and the concave boundary that belongs to other circles: 

d x D(0 ig ) = dB(0,g)nD(0,g), (5) 

d v D(0,e)=dD{Q,e)\d x D(0,Q). (6) 

We note that d x D{ 0, g) = dB{ 0, g)C\V (0). By the Inscribed 
Voronoi Domain Lemma, the vertices of V (0) are all at the 
same distance from 0. This implies that for tl < Q < Rl, 
the convex boundary consists of 2, 4, or 6 circular arcs that 
alternate with the same number of circular arcs in the con¬ 
cave boundary. 



Figure 2: The shaded partial disk D(0, p) with its boundary divided 
into a solid convex portion and a dotted concave portion. 


Angles. Recall that the convex boundary consists of at 
most three pairs of arcs, and let (pi(g) be the angle of each of 
the two arcs in the i-th pair, for i = 1, 2, 3. The total angle 
of the convex boundary is 

3 

<Me) = E 2<pi(e), (7) 

i=l 

and the total angle of the concave boundary is 27r — 'I 1 rfg). 
We have ^l{tl) = 27t and <Pl{Rl) = 0, and between these 
two limits, the function is continuous and mono tonic ally de¬ 
creasing. 

Lemma 4 (Monotonicity) Let L be a lattice in R 2 . Then 
[tl,Rl\ ~t [0, 27t] is continuous, with ^(pi) > 
^l(Q 2 ) whenever pi < p 2 . 

Proof. The continuity of the function follows from the fact 
that dB( 0, g) intersects dV (0) in at most a finite number of 
points. 

To prove monotonicity, we recall that <tv, (p) is the total 
angle of d x D( 0, g) = dB{ 0, g) fl V(0). The Voronoi do¬ 
main is a convex polygon with 0 £ V(Q). Drawing circles 


with radii pi < p 2 centered at 0, we let 0 < 9 < 2n and 
write p\{9) and P2{9) for the points on the circles in direc¬ 
tion 9. Either both points belong to V (0), both points do not 
belong to V(0), or p\(9) £ V(0) but p 2 ($) ^ V(0). The 
fourth combination is not possible, which implies ^(pi) > 
$ L (p 2 ). To prove the strict inequality, we just need to ob¬ 
serve that there is an arc of non-zero length in d x D( 0, pi) 
such that the corresponding arc in dB(0, P 2 ) lies outside 
V (0) and therefore does not belong to d x D( 0, P 2 ). □ 

Area. The probability that a random point belongs to ex¬ 
actly one disk is the area of D( 0, p) over the area of V(0). 
The latter is a constant independent of the radius. We will 
prove shortly that the former is a unimodal function in p with 
a single maximum at the radius p = pl that balances the 
lengths of the two kinds of boundaries; see Figure 3. We call 
Pi the equilibrium radius of L. Write A^: [tl,Rl\ —> R 
for the function that maps p to the area of D( 0, p). 



Figure 3: Increasing the radius grows D( 0, p) along the convex 
boundary and shrinks it along the concave boundary. 

Lemma 5 (Equilibrium Radius) Let L be a lattice in R 2 . 
The function Al : [7'l, Rl] —t R is strictly concave, with a 
unique maximum at the equilibrium radius gL that satisfies 
$ l{QL) = 7T- 

Proof. Recall that 4>i(p) is the total angle of the convex 
boundary of D{ 0, p), and 27r — ^^(p) is the total angle of 
the concave boundary. When we increase the radius, the par¬ 
tial disk grows along the convex boundary and shrinks along 
the concave boundary. Indeed, the derivative is the difference 
between the two lengths: 

^(e) = e[<Me)-(27T-<Mp)] (8) 

= 2p[$ L (p)-7r]. (9) 

The derivative vanishes when ( P/f g) = 7 r, is positive when 
^l(q) > 7r an d negative when ^(p) < tt, as claimed. □ 

4 Optimality of the Regular Hexagonal Grid 

In this section, we present the proof of our main result. After 
writing the probability that a random point lies in exactly one 
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disk as a function of the radius, we distinguish between three 
cases, showing that the maximum is attained at the regular 
hexagonal grid. 

Probability. Given a lattice L in JR 2 , we write < Ql < 
Hi f°r the packing, equilibrium, and covering radii. Recall 
that the probability in question is 

PM = |4tSfe < 1() ) 

in which L = L(a,b) and 7 is the angle between a and b. 
Recall furthermore that the convex boundary consists of at 
most six arcs, two each with angle ip 1 , tp 2 , ¥> 3 , in which we 
set the angle to zero if the arc degenerates to a point or is 
empty. 

Lemma 6 (Equilibrium Area) Let L = L(a, b ) be a lattice 
in R 2 with angle 7 between a and b. Then 

3 

Pl(Ql) = || 0 ||Wsin 7 '5Z Sin ^- ( n ) 

i =1 

Proof. Recall that ||ez|||| 6 || sin 7 is the area of V(0). Let 
Ai n be the area of B{ 0, ql) IT V(0), let A out be the area 
of B( 0, Ql)\V (0), and note that A- ln — A out is the area of 
D{ 0, q l ). Since A in + A out = g 2 L n. we have A in - A out = 
q 2 l n — 2A out . The portion of 5(0, ql) outside the Voronoi 
domain consists of up to three symmetric pairs of disk seg¬ 
ments, with total area 


3 


X-A 2 

Aout = 2 2_^ ^ (tpi - sin (Pi) 

2—1 

(12) 

= ell f -^ sin ^ . 

(13) 


i=l 


in which the second line is obtained using JV =1 (fit = f 
from the Equilibrium Radius Lemma. The probability is 
Aj n — A out divided by the area of the Voronoi domain: 

PM = SSfe. (14) 

Together with (13) this implies the claimed relation. □ 

Case analysis. We focus on the primitive case in which 
the Voronoi domain is a hexagon, considering the non¬ 
primitive case a limit situation in which two of the edges 
shrink to zero length. Let a,b £ R 2 be generators of the 
lattice satisfying the condition in the Non-obtuse Generators 
Lemma, set c = a — b. and assume INI < IHI < ||c||. Re¬ 
call that these three lengths are the distances between parallel 
edges of the hexagon. Further notice that we have Tl = ^ 
for the packing radius and Rl > 77 - for the covering ra¬ 
dius. As before, we write 77 for the angles of the arcs of 
d x D{ 0, q), and we index such that ipi > tp 2 > P 3 - 


Case 1: < Ql < 77 . Then 77 > 0 and <p 2 = 7*3 = 0. 

Case 2: W < ql < Then pi > tp 2 > 0 and tp 3 = 0. 

Case 3: ^ < ql < Rl- Then tpi > > <^3 > 0. 

For example the configuration depicted in Figure 2 falls into 
Case 2. Using the expression for the probability in the Equi¬ 
librium Area Lemma, we determine the maximum for each 
of the three cases. Here we state the results, referring to 
Appendix A for the proofs. By the probability we mean of 
course the probability that a random point belongs to exactly 
one disk. 

Lemma 7 (Two Arcs) In Case 1, the maximum probability 
is attained for ||6|| = -\/2|M|, 7 = arccos ^=, and ql = 

||a||/v/2, which gives Pl{ql) = 0.755 .... 

Lemma 8 (Four Arcs) In Case 2, the maximum probability 
is attained for ||6|| = ||a||, 7 = arccos(v / 2 — 1), and ql = 

ll a ll\/1 — l/\/2> which gives Pl{ql ) = 0.910 .... 

Lemma 9 (Six Arcs) In Case 3, the maximum probability 
is attained for ||6|| = ||o|| = ||c||, 7 = f an d ql = 
||a||/(2 cos jt;), which gives Pl(ql) = 0.928_ 

Note that the lattice in the Six Arcs Lemma is the regular 
hexagonal grid. Comparing the three maximum probabili¬ 
ties, we see that the regular hexagonal grid gives the global 
optimum; see Figures 4. For this lattice, we get Qu such that 



Figure 4: Left: for %/2||a|| = ||fc|| < ||c||, the angle 7 is between 
arccos and Right: for ||a|| = ||b|| < ||c||, the angle 7 is 
between ” and ^ ■ In both cases, the probability increases as the 
angle decreases, attaining its maximum at the minimum angle. 

each disk overlaps with six others and in each case covers 
30° of the bounding circle: the 12-hour clock configuration 
in the plane. This implies the Main Theorem stated in Sec¬ 
tion 1. We further illustrate the result by showing the graph 
of the function that maps ||a||/|| 6 || and the angle 7 to the 
probability at the equilibrium radius; see Figure 5. 

5 Discussion 

The main result of this paper is a proof that the 12-hour clock 
configuration of disks in the plane maximizes the probability 
that a random point lies in exactly one of the disks. Other 
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Figure 5: The graph of the function that maps 0 < IMI/IHI < 
1 and arccos Jpjjj- < 7 < 90° to the probability that a random 
point lies in exactly one disk. The thus defined domain resembles a 
triangle and decomposes into three regions corresponding to Cases 
1, 2, 3. The regular hexagonal grid is located at the lower left comer 
of the domain. 

criteria favoring configurations between packing and cover¬ 
ing can be formulated, see [7], and it would be interesting 
to decide which one fits the biological data about DNA or¬ 
ganization within the nucleus best. There are also concrete 
mathematical questions related to the work in this paper: 

• Is the 12-hour clock configuration optimal among all 
configurations of congruent disks in the plane? 

• What is the optimal lattice configuration of balls in R 3 ? 

To appreciate the difficulty of the second question, we note 
that the FCC lattice gives the densest packing [6], while the 
BCC lattice gives the sparsest covering [1], Does one of 
them also maximize the probability that a random point lies 
inside exactly one ball? 
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A Proofs 

In this appendix, we give detailed proofs of the three Arc 
Lemmas. As described in Section 4, the three lemmas add 
up to a proof of the Main Theorem stated in Section 1 of 
this paper. We begin with a few relations that will be useful 
in all three proofs. Given a triangle with edges of lengths 
||a||, || 6 ||, ||c|| and angle 7 opposite the edge c, the law of 
cosines implies 

\\c\\ 2 =\\a\\ 2 + \\b\\ 2 -2\\a\\\\b\\cos r (15) 

Assuming ||a|| < || 6 || < ||c|j, the angle 7 is at least as large 
as each of the other two angles. From (15) together with 
||fr|| < ||c||, we get COS7 < Jpjy- As justified by the Non- 
obtuse Generators Lemma, we may assume the triangle is 
non-obtuse, which implies 


Proof of the Four Arcs Lemma. Case 2 is defined by 

< ql < 7 p, which implies (pi > > 0 and = 0. 

The probability at the equilibrium radius is therefore 

P l{Ql) = | [a| | 2 b fsin 7 ' ( sin V’t + sin^ 2 ). (23) 

To get a handle on the maximum of this function, we first 
write the the sum of sin ipi and sin ip 2 and second the equi¬ 
librium radius in terms of other parameters. Using cos 2 a = 
cos 2 a — sin 2 a and cos 2 a + sin 2 a = 1, we get cos p 2 = 
2 cos 2 iy- — 1, and since tpi + ip 2 = §, we have sinyq = 
cos p 2 . Recalling (18), we get sin <p 1 = \\b\\ 2 / {2q 2 l ) — 1, and 
recalling (17), we get siny> 2 = ||a||“/(2p 2 ) — 1. Adding the 
two relations gives 

, • ||a|| 2 + ||6|| 2 o rsA\ 

sin (fix + sm <p 2 = o 2 - 2 . (24) 


arccos lp{f < 7 < f • (16) 

The range of possible angles is largest for ||a|| = ||b|| where 
we get 60° < 7 < 90°. Furthermore, we write the angles 
ifi of the arcs in the convex boundary of the partial disk in 


terms of the edge lengths and the radius: 


cos £1 _ M 

(17) 

cos £2 _ M 

CUb 2 — 2q ’ 

(18) 

v?3 _ INI 

2 2 Q ’ 

(19) 


The first relation holds provided 77 < g < R/ , and sim¬ 
ilar for the second and third relations. Finally, we note that 
scaling does not affect the density of a configuration. We can 
therefore set || 6 || = 1 , which we will do to simplify compu¬ 
tations. 


Proof of the Two Arcs Lemma. Case 1 is defined by 

7 P < ql < which implies > 0 and (p 2 = p 3 = 0 . 
Since sin p 2 — sin (p$ = 0, the probability at the equilibrium 
radius simplifies to 


Pl{ql) — || 0 |||| b j| L s i n7 • sin ipi 
= M 

||&||sin 7 ’ 


( 20 ) 

( 21 ) 


where we get the second line by combining pi = | with 
(17) to imply ql = For the remainder of this proof, 

we normalize by setting || 6 || = 1. To maximize the proba¬ 
bility, we choose ||a|| as large as possible and 7 as small as 
possible. From ||a.|| /= Ql < we get ||a|| < l/x/2, 
and from || 6 || < ||c|| we get 7 > arccos The two pa¬ 
rameters can be optimized simultaneously, which gives 


Pl(ql) 



0.755.... (22) 


To find a substitution for the equilibrium radius, we begin 
with (18), use ff = j — and finally apply cos(cc + /?) = 
cos a cos /? — sin a sin /3: 


M 

2 ql 


= COS 


a) 


1 

V2 


(cos ^ + sin . 


(25) 

(26) 


Next, we substitute the two trigonometric functions using 
(17) and sin 2 a = 1 — cos 2 a. Simplifying the resulting 
relation and squaring it, we get 


Q 2 L = h{\M 2 + m 2 -V2\\a\\\\b\\). (27) 

Plugging (24) and (27) into the equation for the probability 
and normalizing by setting || 6 || = 1 , we get 


Pl(ql) 


2U2||q|H|a|| 2 -l 
IIall si«7 

2 j n/2||q|| — 11 a 11 2 1 


1- v^- 


klr+i 


(28) 

(29) 


where we maximize to get the second line by choosing 7 as 
small as possible. Specifically, 7 is implicitly restricted by 
W- < ql < 77 , so we can use Aq 2 l < ||c || 2 together with 
(15) and (27) to get 

COS'Y < yp) — + (30) 

cos 7 ^ V ^ 2]]q|| \\b\\ ■ 

Checking with the Maple software [11], we find that the 
right-hand-side of (29) increases in [0,1] attaining its max¬ 
imum at ||a|| = 1. We therefore get 7 = arccos (\/2 — l) 
from (30), p 2 = 1 — l/y/2 from (27), and 


Pl{ql) 



0.910... 


(31) 


from (29). 
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Proof of the Six Arcs Lemma. Case 3 is defined by 

< Ql < Rl , which implies <£>i > p >2 > <£>3 > 0. Start¬ 
ing with the expression for the probability given in the Equi¬ 
librium Area Lemma, we first express the sin </?,; in terms of 
the other parameters: 


Sint/?! = 




2gf, 


_ l|b|lv / 4^HI 


sin <£> 2 = 


sin <£>3 = 1 — 


2 ei 


\\4 V^PlW+WH v / 4elHHi 


set 


(32) 

(33) 

(34) 


To get (32), we use sin 2a = 2 sin a cos a with a = 
together with (17). To get (33), we use the same trigonomet¬ 
ric identity with a = , together with (18). To get (34), 

we use the equilibrium condition together with sin £>3 = 

sin(| - <£>1 - <£> 2 ) = cos((£>i + ip 2 ) = 1 - 2 sin 2 yi + y2 , 
and finally substitute sin(a + (3) = sin a cos f3 + sin /? cos a, 
with a = and (3 = ^. 

To do the same for sin 7 , we take the cosine of both sides 
of the equilibrium condition, which is ^ = j — ^ ^. 

Writing <7 = cos ^ and s,; = sin for * = 1,2,3, and 
applying standard trigonometric identities, we get 

° 3 = ^73 [ ClC2 - SlS 2 + Sl c 2 + c lS 2 }- ( 35 ) 

c\ = \ + ( 2 ciC 2 - cO^l - cf 

+ (2c 2 c 2 -c 2 )-y/l-c 2 . (36) 


Using (17), (18), (19) and substituting ||c || 2 using (15), we 
get the following relation after a few rearrangements: 
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(37) 


Using cos 2 7=1 — sin 2 7 , we can substitute sin 7 in the for¬ 
mula for Pl(ql)- We thus arrived at a relation that gives the 
probability in terms of |ja||, || 6 ||, and ql only. While being 
lengthy, this relation is readily obtained by plugging (32), 
(33), (34), and (37) into (11). We therefore take the liberty 
to omit the formula here and refer the interested reader to the 
website of the second author of this paper 1 . 

It remains to determine the parameters that maximize the 
probability. To simplify this task, we normalize by setting 
ll&H = 1. The probability is thus a function of two variables, 
||a|| and ql. Using the Maple software, we compute the two 
partial derivatives, 9P/,/9||a|| and OPlI^Ql- Setting both 


*A Maple file with the main steps in the formulas related with 
this paper is available at http: //mabelih9 . wix. com/mabelhome# 
!publications/cee5. 


to zero, we get ||a|| = 1 matched up with three radius values: 

ql = | (/6- V^) = 0.517..., (38) 

ql = -i + C7-3 = °.5 82 ..., (39) 

q l = 4 (\/6 + \/2) = 1.931..., (40) 

with C = 3 + 2\/2. Setting j|a|| = ||6|| = 1, the covering 
radius depends only on the angle 7, which ranges from | to 
|. It is largest for 7 = |, where R L = \/2/2 = 0.707.... 
The radius we get at the third root is larger than that and can 
therefore be excluded. 



Figure 6 : The quadrangle in the plane defined by the length of a 
and the radius is shaded. Along its boundary, we encounter three 
local minima (two corners and a point along the right edge) and 
three local maxima (a corner and a point each along the lower edge 
and the right edge). 


The maximum probability is attained at one of the two re¬ 
maining roots or along the boundary of the domain region 
that corresponds to Case 3. As illustrated in Figure 6 , we 
simplify the computation by taking a quadrangle that con¬ 
tains this region. The quadrangle is defined by 
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(41) 

(42) 


in which the first interval follows from the bounds that define 
Case 3, and the second interval is obtained by limiting the 
radius by the covering radius of the configuration in which 
a and b enclose its maximum angle. The maximum angle 
for Case 3, 7 max* corresponds to the minimal angle for Case 
2 derived from (30). The upper boundary on g is a convex 
curve and so we replace it with the straight line connecting 
its extremes. 
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We evaluate the probability at the four vertices and, us¬ 
ing the Maple software, at the roots of the derivatives along 
the four edges. As shown in Figure 6, we find three lo¬ 
cal minima alternating with three local maxima along the 
boundary of the quadrangle. The maximum and minimum in 
the interior of the right edge coincide with the roots in (38) 
and (39). Among the three local maxima, the probability is 
largest at (38), which is characterized by ||a|| = ||b|| = 1 and 
Ql = \ (\/6 — V2) = 1/(2 cos ^). This gives P l (ql) = 
0.928 ..., as claimed in the Six Arcs Lemma. Indeed, plug¬ 
ging the values into the equilibrium condition of Case 3, we 
get ||c|| = 1 , which shows that the probability is maximized 
by the regular hexagonal grid. 




